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1. Introduction 

The simulation of final states in the collision of highly energetic particles at current 
and future colliders is largely based upon parton shower evolution models. The 
theoretical foundation of these models is the interpretation of Sudakov form factors 
as no-emission probabilities 0. For hadronic collisions, the space-like evolution of 
initial state partons from the hard scale Qh of the partonic subprocess and relatively 
small X is 'guided' by the ratio of parton distribution functions (pdf's) [If. In this 
paper we would like to address the effect of uncertainties in the Sudakov form factors 
as a result of pdf uncertainties 0, Q . We shall compare this effect to the uncertainty 
due to the choice of scale in the strong coupling constant as- We choose to use 
the Sudakov form factors that are going to be used for the initial state evolution 
in future versions of Herwig++ 0. The conclusion should qualitatively remain the 
same for the Sudakov form factors that results from different kinematical variables in 
the programs HERWIG and PYTHIA[|| or the new program SHERPA ||. They 
may, however, differ in details as the different kinematics will lead to emphasis on 
different regions of evolution phase space. 

Furthermore we limit our discussion to the case of spacelike evolution. The 
estimated size of uncertainties due to effects like scale variation in as{Q) should 
be similar in the case of timelike evolution. As the pdf's only appear in spacelike 
evolution we can study all effects of our interest here. Another limitation is the 
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discussion of the uncertainties due to terms in the Sudakov form factor that are 
formally sub leading in the sense of the next-to-leading (soft and/or collinear) log- 
arithmic approximation. Formally we may add or subtract terms in the splitting 
functions that give rise to only non-logarithmic terms after integration over q (see 
below). Adjusting these terms can be crucial for matching of matrix elements and 
parton showers in the CKKW scheme |TD|, |TT| as discussed in [TB|. The effect, 



however, appears to be mostly of technical importance close to the matching scale 
inherent to this kind of matching scheme. 



2. Sudakov form factor for space— like branchings 

The Sudakov form factor for spacelike backward evolution of a parton a from the 
hard scale gmax down to some scale q can be written as p|, 0] 



Q'maxi -^j Q'O / 



exp 



.x,qo} 



(2.1) 



The sum on the right hand side (rhs) is over all possible splittings into partons of 
type b and 
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Figure 1: 

branching b 



Kinematics for space-like 
ac. 



We assume a resolution scale go below 
which the evolution terminates and further 
splittings would remain unresolvable. The 
limits of the ^-integration, zq and Zi de- 
pend imphcitly on q and hence the go 
dependence^ of 2ba and Sa- x is the light 
cone momentum fraction of the parent par- 
ton with respect to the originating hadron. 
This value is initially selected by the hard 
subprocess. x' = x/z is the light cone mo- 
mentum fraction of the new parent after the first space like branching and so forth. 
Pbaiz^q^) is the unregularized collinear splitting function, which, in the case of the 
evolution of massive partons, may also depend on the branching scale (f. Note, that 
the splitting function is regularized as we use explicit cutoffs in the phase space for 
soft gluon emission. /a(a;, g^) is the parton distribution function (pdf) of a parton of 
type a inside a proton^. 

The integral in Eq. ( p.2|) is written in a fairly symbolic way and holds for different 
models. We shall rewrite it for the specific case of the kinematic variables that 

^Somc authors prefer to stress this circumstance by writing Sha as a ratio of two exponentials 
with explicit go dependence in the g-integration. The two notations are, of course, equivalent. 

^The conclusions will remain unchanged when considering antiprotons. The case of pious or 
resolved photons in the initial state is beyond the scope of this study. 
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are used in Herwig++ |14|. There, the evolution variables z,q are interpreted in a 
Sudakov basis. Consider the branching in Fig. |l|. Each momentum is written as 



Qi = ttip + PiU + q±i 



(2.3) 



where p and n are suitable forward and backward light like vectors, typically chosen 
along the axis of the original hadrons. Then we define our space-like branching 
kinematics via 

ai = , q±i = . (2.4) 



Z Zi 

Here, the relative transverse momentum p±i is given in terms of the evolution vari- 
ables q and z as 



Pi. 



(1 



\2~2 



[2.5) 



Qg is the cutoff parameter of the parton shower. Choosing the argument of as{Q) 
a.s Q = {1 — Zi)qi we may now rewrite the integral ( p.2D as 



ymax 



as[{l - z)q] x'fb{x',q^ 

n2\ 



-3r / dz 



27r 



xfa{x,q 



(2.6) 



The implicit introduction of 
the parton shower cutoff Qg in 
( ^.5| ) is reflected in the step func- 
tion G(P.S.). Hence we simply 
extend the integration region in 
z from to 1 and apply this cut 
on the available parton shower 
phase space (see also Fig. Q) . How- 
ever, the resolution scale go of the 
parton shower has disappeared and 
was implicitly replaced by Qg. Now 
every given initial condition gmax, x 
may result in different values of 

the resolution scale in terms of Figure 2: Available phase space for space-like 
a smallest available scale. The branching, 
shape of the available phase space 

is, of course, still determined by Qg in some intuitive way such that small cutoffs lead 
to a larger available phase space and a lower resolution scale. We set Qg ^ 1 GeV if 
not stated otherwise^. 

In practice, the available phase space is determined by some more factors. As 
we require a real transverse momentum we find the limits of the z integration for a 




^More precisely we choose the HerwigH 
to but shghtly smaUer than 1 GeV. 



parameter 6 — 2.3 GeV (cf. 0). This leads to Qg close 



3 




given q. As x' should never exceed 1, x itself is the lower limit of z. Therefore, 

(2.7) 

Depending on the choice of how to evaluate as{Q) when Q becomes very small 
we have an additional constraint on the phase space. Let Qo be the scale where 
non-perturbative effects become important. Typically Qq is somewhat larger than 
but of the order of Aqcd- In order to model non-perturbative effects, we may simply 
freeze the value of as, O-siQ < Qq) = C(s{Qo)- We set Qq = 0.75 GeV and freeze as 
by default. This gives us a behaviour of very similar to the one in the MRST2001 



parametrization as we find it in |T^. If we put as{Q < Qo) = we have an 
additional limitation of phase space. In Fig. ^ we show as an illustrative example 
the available phase space for q ^ qg spacelike branching with gmax = 4 GeV and 
X = 0.02 with its limiting regions. Lines of constant as for Qq/2, Qq,2Qq are also 
shown. 

We apply an extra cut at q < Qg in the rare cases where there is still phase space 
available. As the pdf 's become unstable for Q values close to the non-perturbative 
domain we choose to freeze the pdf 's at Qfreeze = 2.5 GeV. The choice of this param- 
eter is based on experience with HERWIG^. 

In the following we consider the Sudakov form factor for specific branching types 
6 — >■ ac only. Formally we rewrite ( p.l|) as 

Qg) = YlSbaiq, Qg) (2.8) 

b 

where Sha = exp —Iba- Then we only consider specific Sba- For the parton shower 
evolution the Sudakov form factor is interpreted as the probability that the parton a 
remains unresolved upon evolution from q^ax down to q when the minimal resolution 
is go- With the probability 1 — Sba that there is any branching between gmax and q 
we find the branching probability density 

Vbaiq, gmax; X, Qg) = ll^{q; X, Qg)Sba{q, Qg) , (2.9) 

where X'f^^{q;x,Qg) is the integrand of ( |2.6| ) with respect to the q integration. The 
scales of the next branching q are distributed according to Vbaiq, Q'max; x, Qg). This is 
the quantity we would like to consider in order to estimate uncertainties of different 
types as it is directly implemented in parton shower Monte Carlo programs. We 
should note that Vbaiq, qmax] x, Qg) is normalized to 

/ Vbaiq, O'max; X, Qg) dq = 1 - Sbaiqo, ^max; X, Qg) . (2.10) 

This is typically different from 1 as we have a chance Sbaiqo, qma.x',x,Qg) to find no 
branching at all. Their sum is obviously one. 

^In HERWIG the pdf 's are regularized by the parameter QSPAC with the defauh value 2.5 GeV. 



-4- 



3. Numerical study 



For numerical purposes we map the integration variables q and z onto new variables 
t, y as 



t = \nq\ y = \n^, (3.1) 
in order to absorb the leading poles into the integration variables, 



dt = ^, dy=^^^. (3.2) 
Z{1 — z) 



We carry out the integrations in ^^^(g; x, Qg) and Iba{(j, Qmax] x, Qg) numerically with 
adaptive Gaussian or Monte Carlo algorithms |T6|, [1^. Different packages have 
been used for consistency checks. Particular care has been taken to obtain numerical 
results with a typical relative error of at least one order of magnitude below the 
resulting pdf errors as the former ones would mask the latter ones otherwise. 

In the following we address the uncertainties of Sudakov form factors that arise 
from various sources. As we do not aim at a study of Sudakov form factors beyond 
leading logarithmic accuracy, the sources of uncertainties are: 

1. the ratio of pdf's, which may even be pdf's of different flavour in the case of 
q ^ gq and g ^ qq branching, 

2. parametrization of the running of 0:5 may play an important role, 

3. modelling of as in the non-perturbative region below our choice of Qo- 

Recently, the uncertainties of parton distribution functions have been considered for 
various observables and the question has been raised how much of an influence they 



would have on Monte Carlo showering algorithms ||T8[| . 

In order to assess the relevance of uncertainties arising from the parton distribu- 
tion functions we basically follow the procedure described in and We obtain 
the integral once for the central pdf fit and once for all the error members. The 
variations from each pair of error pdf's are added in quadrature and thus we obtain 
an error estimate 6Iha for the integral from which we can derive the uncertainty in 
the probability distribution |2.9|. All calculations have been performed with the pdf 



error distributions from MRST and CTEQ P, H as they are provided in |TP 



The parametrization of the running of as{Q) may give important differences in 
the result. Particularly as different orders of the running result in large deviations at 
smaller values of Q. In Fig. |^ we show our default parametrization with and without 
freezing as well as the MRST parametrization of as{Q)- By default we use the two- 
loop as that is implemented in Herwig++. We choose to freeze as at Qo = 0.75 GeV 
as discussed above. This parametrization is numerically very close to that of MRST. 
The parametrization of the 05 in the CTEQ distribution is very similar but does not 
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model as at small scales Q in some particular way. We therefore do not check this 
parametrization explicitly. For MRST we compared results from a computation with 
the given pdf set where by default we use our own parametrization of as with the 
same computation where the ^5 is taken from the distribution of the pdf set itself. 

As the results seemed to be fairly 
sensitive on the choice of a^-para- 
metrization we have also checked the 
effect of setting as to zero below Qq. 
This may have drastic consequences 
as one can see with the help of Fig. |[ 
There, we plotted lines of constant 
Qo once for the central value and once 
the values resulting from rescaling by 
factors of 2. We see that we may cut 
out regions in phase space that actu- 
ally give very large contributions to 
the integral Iba- For estimates of this 
scale uncertainty we compute the Su- 

dakov form factors and the branching probability densities with variation of the scale 
by a factor two. 
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Figure 3: Different parametrizations as{Q)- 



4. Results 

We consider the Sudakov form factors and branching probability densities in the 
following. We choose different initial conditions qmax,x for the evolution and focus 
on a particular type of branching at a time. The plots contain three panels and show 
Sfya on top, Vba in the middle and some error information on the bottom panel. 

4.1 Different initial conditions 

We first consider different initial conditions and branching types. In this case we 
always use MRST partons. We plot Sudakov from factor and branching probability 
density for a number of different initial conditions. The error box in the lowest 
panel always compares the error from the numerical integration with the total error 
that also contains the pdf error itself. This is to control that the numerical error is 
negligible compared to the pdf error. This is always the case as the numerical error 
is typically two orders of magnitude smaller than the pdf error. 

In Fig. ^ we show the uncertainties in the case oi q ^ qg spacelike branching. 
We pick d quarks in particular. Clearly, the estimated pdf error is by far smaller than 
the one from varying the as scale by a factor of two. Generally, the pdf uncertainty 
is largest at small initial scales gmax = 10 GeV and somewhat larger at small x. This 
goes along with the fact that the quark densities are known better at large x than at 
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small X. We do not show any plots at intermediate values of x and q as they don not 
exhibit any new features. The plots shown in Fig. ^ are those with a large visible 
effect. 

For g ^ gg branchings we get a slightly different picture in Fig. |^. The pdf 
uncertainties are much larger in this the gluon density is not as much con- 

strained as the valence quark density. This is particularly true at large x. Going to 
small X = 10~^ at (jmax = 10 GeV the pdf uncertainty is very large (the largest we 
encounter in this study) and sometimes even larger than the uncertainty from as 
scale variation. This is still the case if one goes to a larger initial q^a.^ = 30 GeV but 
the effect is still largest at small target values of q. For even larger gmax (not shown) 
we observe the same trend. 

We also consider the more unlikely branchings g gq and g ^ qq m Fig. ^ 
The Sudakov form factor for these branchings is typically close to one. The pdf 
uncertainties can be sizable as well, particularly when we involve a gluon at small 
values of x where we also had the largest uncertainty in the case of g ^ gg splitting. 
Still, we have selected situations with particularly large pdf uncertainties. Generally, 
we have the larger error from as scale variation. 

In Fig. 1^ we replace the initial valence quark a = d with a sea quark a = d. 
Now we can directly observe and compare the sea-like behaviour of the d as opposed 
to the valence d-quark. Especially at larger a; = 0.1 the evolution clearly prefers 
valence quarks (note, that the branching probability density is normalized to the 
total probability of such a branching to occur). We can directly compare to the 
picture at small x = 10~^ where the two distributions begin to coincide as the 
valence structure loses its dominance. The opposite happens of course for evolution 
away from the d/d into a gluon. The evolution of the d is clearly suppressed. As 
far as the errors due to pdf 's are concerned we observe that the uncertainty of the 
sea quark density at large x becomes visible as it does for the gluon. This is highly 
emphasised at the low end of the evolution. 

4.2 MRST vs CTEQ 

As an obvious cross check we repeated all computations also with the CTEQ6m error 
set. In Fig. |^ we show some of the previous plots again, without as uncertainties, 
comparing the error bands of the two parton distribution packages. We find that 
the distributions are compatible but the error estimates from CTEQ are smaller, 
particularly for gluons at small x. Again, we have picked cases that show clear 
differences. Generally, both packages tend to give similar (and compatible) results. 

4.3 Parametrization ot as 

As already pointed out in the discussion of the available phase space, the effect of 
different parametrizations of as{Q), especially at low Q can have important effects. 
In Fig. ^ we show the effect of using the as parametrization that is provided via 
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the distribution of MRST ifself with our parametrization. We have chosen a very 
similar regularisation for non-perturbative values of as as MRST have. Therefore 
the results are very similar, nevertheless they are sinsitive to the small differences in 
the parametrization of as- However, the effect of using a parametrization that sets 
as to zero below e.g. Qo = 0.75 GeV leads to a very different result as effectively a 
part of phase space is modified in which X^)a picks up very large contributions. This 
is of particular importance when we use Q/2 as this cuts out a lot of phase space. 
In the other two cases the results closely follow those of our default parametrization 
until they hit the phase space boundary. 



4.4 Larger cutoff 

As a final subject we consider the spacelike evolution of light quarks or gluons at 
large x = 0.2 and fairly high initial scale gmax = 250 GeV. This situation can give 
us some information about the initial state radiation in top production at hadron 
colliders, ti pairs are predominantly produced from light quark scattering at large x. 
The relevant initial scale gmax is given by a typical invariant mass in the t-channel as 
the colour connection is between one initial light quark and one t-quark [Q. This 
gives us a typical initial scale gmax = 250 GeV. 



In Fig. |T0| we have plotted the relevant Sudakov form factors and branching 
probability densities for such a situation. We have increased the cutoff scale to 
higher values Qg = 2.5 GeV and/or 10 GeV. This allows us to directly estimate the 
effect of the radiation of an extra gluon at a fairly large (visible) scale. We find that 
the pdf uncertainties remain fairly low in the Sudakov form factors themselves. The 
size of the branching probability distribution is affected to some extent but the shape 
is hardly modified. Again, the uncertainty from the variation of as is much larger. 



5. Conclusion 



We have studied the effect of various sources of uncertainties on the parton shower 
evolution. We therefore particularly considered the branching probability density as 
this is the quantity that drives the evolution. We generally find that the effect of 
pdf uncertainties is small compared to uncertainties in the QCD description of the 
evolution itself that we estimate via scale variation in the strong coupling. There 
are, however, particular regions in phase space where the pdf uncertainties are very 
large. This is the case for the evolution from very small values of x and also at large x 
when sea quarks or gluons are involved. Considering the initial conditions in q we find 
that the pdf uncertainties are generally small when the initial conditions are in the 
large q regime. In the course of the parton shower evolution, however, all showers 
have to evolve through the small q in some way. Besides, we emphasise that the 
evolution strongly depends on the modelling of as in the transition region between 
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perturbative and non-perturbative QCD. This uncertainty is, however, usually taken 
as a sensitivity that is exploited in tuning Monte Carlo showers to experimental data. 

In summary, we find that one should be aware of pdf uncertainties in the evo- 
lution when one is confined to a phase space of fairly large x and relatively small 
initial q. 

Acknowledgments 

I would like to thank J. Huston for fruitful conversations and for raising this topic 
at the HERA/LHC workshop. I also thank Peter Richardson, Mike Seymour and 
Bryan Webber for valuable comments. 

References 

[1] G. Marchesini and B. R. Webber, Nucl. Phys. B 310, 461 (1988). 

[2] T. Sjostrand, Phys. Lett. B 157, 321 (1985). 

[3] R. K. Ellis, W. J. Stirling and B. R. Webber, "QCD and Collider Physics," 
Cambridge Monogr. Part. Phys. Nucl. Phys. Cosmol. 8 (1996) 1. 

[4] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 28 

(2003) 455 [arXiv:hep-ph/0211080]. 

[5] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 
0207 (2002) 012 [arXiv:hep-ph/0201195]. 

[6] S. Gieseke, A. Ribon, M. H. Seymour, P. Stephens and B. Webber, JHEP 0402, 005 

(2004) [arXiv:hep-ph/0311208]. 

[7] G. Corcella et al, "HERWIG 6.5 release note," arXiv:hep-ph/0210213; 

G. Corcella et al, "HERWIG 6: An event generator for hadron emission reactions 
with interfering gluons (including supersymmetric processes)," 
JHEP 0101 (2001) 010 [arXiv:hep-ph/0011363]. 

[8] T. Sjostrand, L. Lonnblad, S. Mrenna and P. Skands, "PYTHIA 6.3 physics and 
manual," arXiv:hep-ph/0308153; T. Sjostrand, P. Eden, C. Priberg, L. Lonnblad, 
G. Miu, S. Mrenna and E. Norrbin, "High-energy-physics event generation with 
PYTHIA 6.1," Comput. Phys. Commun. 135 (2001) 238 [arXiv:hep-ph/0010017]. 

[9] T. Gleisberg, S. Hoche, F. Krauss, A. Schalicke, S. Schumann and J. C. Winter, 
JHEP 0402 (2004) 056 [arXiv:hep-ph/0311263]. 
http : //www . physik . tu-dresden . de/~krauss/hep/. 

[10] S. Catani, F. Krauss, R. Kuhn and B. R. Webber, JHEP 0111 (2001) 063 
[arXiv:hep-ph/0109231]; 



-9- 



[11] F. Krauss, JHEP 0208 (2002) 015 [arXiv:hep-ph/0205283]. 

[12] S. Mrenna and P. Richardson, JHEP 0405 (2004) 040 [arXiv:hep-ph/0312274]. 

[13] F. Krauss, A. Schalicke, S. Schumann and G. Soff, Phys. Rev. D 70 (2004) 114009 
[arXiv:hep-ph/0409106]. 

[14] S. Gieseke, P. Stephens and B. Webber, JHEP 0312, 045 (2003) 
[arXiv:hcp-ph/0310083]. 

[15] W.H. Press, B.P. Flannery, S.A. Teukolsky, W.T. Vetterhng, Numerical Recipes in 
C, Cambridge University Press, 1992. 

[16] http : //www . gnu . org/ software/ gsl/ 

[17] T. Hahn, arXiv:hep-ph/0404043. 

[18] J. Huston, private communication. 

[19] http : / / durpdg . dur . ac . uk/lhapdf / 



- 10 - 




Hcrwig++ spacclikc q ^ qg 



Hcrwig++ spacclikc q ^ qg 




20 



100 GeV, x = 0.1 

H \ 

MRST2001E 

as{Q/2) 

as(20) 




+ 

int+pdf 




40 60 
(j/GcV 



80 



100 




40 60 
q/GcY 



100 



Figure 4: pdf and as uncertainties in spacelike q ^ qg branchings for different initial 
conditions qmax and x. Each plot shows Sha{q) (top), Vba{q) (middle) and the relative error 
SVha/'Pba from numerical integration and pdf uncertainties (bottom panel). The central 
value is shown as solid line, central it pdf error is dot-dashed and as uncertainties from 
rescaling the argument by a factor 1/2 and 2 are shown as dashed/dotted line, respectively. 
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Q'max values. See caption of Fig. § for labelling. 
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Figure 6: pdf and as uncertainties of the branchings q 
x-values. See caption of Fig. ^ for labehing. 
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Figure 7: Comparison of uncertainties in the evolution of valance (d) quarks and sea 
{d) quarks. The lowest panels show the total relative error of the branching probability 
densities. The solid lines are partly identical to the pdf error bands in Figs. Q and |6|. See 
caption of Fig. ^ for labelling. 
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Figure 8: Comparison of pdf uncertainties from MRST2001E (solid) and CTEQ6mE 
(dashed) parametrizations where differences are most significant. Upper and middle panel 
are similar to Fig. ^ The bottom panel shows the relative pdf error of the branching 
probability density. See caption of Fig. Q for labelling. 



- 15 - 



Hcrwig++ spacclikc q qg 



Hcrwig++ spacclikc q ^ gq 




0123456789 10 0123456789 10 

q/GcV q/GcY 

Figure 9: Effect of different as parametrizations in spacelike q ^ qg and q ^ gq branch- 
ings. The soUd line in both plots denotes the previous results, obtained with MRST partons 
and our default parametrization of 05, including the effect of varying the scale by a factor 
of 2. The dashed and dotted lines are results obtained with different parametrizations of 
as- 
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Figure 10: Increasing the cutoff Qg in the q ^ qg and g 
at large cjmax and x. See caption of Fig. for labelhng. 
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